Franceschini 



Simultaneous estimation of attenuation and structure 
parameters of aggregated red blood cells from backscatter 

measurements 



Emilie Franceschini, Francois T. H. Yu, Guy Cloutier 

Laboratory of Biorheology and Medical Ultrasonics 
University of Montreal Hospital Research Center 
2099 Alexandre de Seve (room Y-1619), 
Montreal, Quebec, H2L 2W5, Canada 
Emilie . Franceschini @ crchum . qc .ca 
Guy.Cloutier@umontreal.ca 



Page 1 



Franceschini 



Abstract 

The analysis of the ultrasonic frequency-dependent backscatter coefficient of aggregating red 
blood cells reveals information about blood structural properties. The difficulty to apply this tech- 
nique in vivo is due to the frequency-dependent attenuation caused by intervening tissue layers that 
distorts the spectral content of backscattering properties from blood microstructures. An optimiza- 
tion method is proposed to simultaneously estimate tissue attenuation and blood structure factor. 
In an in vitro experiment, the method gave satisfactory estimates with relative errors below 22% 
for attenuations between 0.101 and 0.317 dB/cm/MHz, signal-to-noise ratios>28 dB and kR <2.7 
(k being the wave number and R the aggregate radius). 

PACS numbers: 43.80.Qf, 43.80.Cs, 43.35.Bf, 43.35.Yb 
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1. Introduction 

Ultrasonic (US) backscattered echoes from blood contain frequency-dependent information 
that can be used to obtain quantitative parameters reflecting the aggregation state of red blood 
cells (RBCs). Recently, two parameters describing RBC aggregation, the packing factor and mean 
fractal aggregate diameter, were extracted from the Structure Factor Size Estimator (SFSE). 1 The 
SFSE is a second-order data reduction model based on the structure factor and adapted to a dense 
medium such as blood. This approach is based on the analysis of the backscattered power spec- 
trum that contains information about the size, spatial organization, concentration and mechanical 
properties of scatterers (i.e. RBCs). The difficulty to apply the SFSE in vivo is that the spectral 
content of backscattered echoes is also affected by attenuation caused by intervening tissue lay- 
ers between the probe and the blood flow. More generally, ultrasound scatterer size estimation 
techniques for tissue characterization (such as liver, kidney, prostate or breast) are facing similar 
challenges. 2,3 To evaluate correctly microstructural parameters, it is thus of major interest to take 
into account tissue attenuation effects. A few groups 2 ' 4 ' 5 developed measurement techniques to 
evaluate the frequency-dependent attenuation in order to compensate a posteriori the backscat- 
tered power spectrum. The goal of this letter is to further develop this strategy for in vivo measures 
of RBC scatterer sizes. We propose to determine simultaneously blood structural parameters and 
total attenuation by using an optimization method, termed the Structure Factor Size and Attenua- 
tion Estimator (SFSAE). 

This method consists in fitting the spectrum of the backscattered radio-frequency (rf) echoes 
from blood to an estimated spectrum by a modified SFSE model. This approach is similar to that 
presented by Bigelow et al. 3 , who estimated the effective radius of tissue microstructure and total 
attenuation from simulated backscattered signals. Herein, in vitro experimental evaluation of the 
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SFSAE is performed. Porcine RBCs were sheared in a Couette flow system, and ultrasonic rf 
echoes were obtained using a 25 MHz center-frequency transducer. Since skin is one of the most 
attenuating tissue layers during in vivo scanning, three skin-mimicking phantoms with different 
attenuation coefficients were successively introduced between the transducer and the blood flow. 
This study shows the SFSAE ability to evaluate three parameters (the packing factor, mean fractal 
aggregate diameter and total attenuation). 

2. Structure Factor Size and Attenuation Estimator 

Blood can be seen as a very dense suspension of red cells. These RBCs cannot be treated as 
independent scatterers since particle interactions (collision, attraction, deformation, flow depen- 
dent motions) are strong. The theoretical model of ultrasound backscattering by blood that we 
developped 1 is based on the particle approach 6,7 , which consists of summing contributions from 
individual RBCs and modeling the RBC interaction by a particle pair-correlation function. Assum- 
ing that the Born approximation is valid (weak scattering), the model proposed in Ref. 1 can be 
modified to predict the theoretical backscatter coefficient from blood: 

BSC theor (k) = ma b (k)S(k)A(k) (1) 

where k is the wave vector, m is the number density of RBCs in blood estimated by measuring 
the hematocrit H by microcentrifugation (m = H/V s , where V s is the volume of a RBC), cr& is 
the backscattering cross section of a single RBC, S is the structure factor describing the spatial 
organization of RBCs, and A is the frequency-dependent attenuation function. The backscatter- 
ing cross-section a b of a weak scattering particle small compared to the wavelength (Rayleigh 
scatterer) can be determined analytically as follows: cr^/c) = l/(47r 2 )/c V/7 2 , where j z is the 
variation of impedance between the RBC and its suspending medium (i.e. the plasma). The struc- 
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ture factor S is by definition the Fourier transform of the pair-correlation function 7 g and can be 
approximated by its second-order Taylor expansion 1 in k as 

S(k) = 1 + m J (g(r) - 1) e~ 2]kr dr « W - ^{kRf. (2) 

In this expression, g(r) represents the probability of finding two particles separated by a distance 
r. W is the low-frequency limit of the structure factor (S(k) |fe_>o) called the packing factor. 7 ' 8 R is 
the radius of 3D RBC aggregates assumed to be isotropic. We introduce D = R/a as the isotropic 
diameter of an aggregate (expressed in number of RBCs) with a the radius of one RBC sphere- 
shaped model of volume V s . The attenuation function A is given by: A(k) = e~ Aa "f, where / is the 
frequency and «o is the attenuation coefficient (in dB/MHz) defined by: a = onei, where ctij 

i 

and ti are respectively the intervening tissue layer attenuations (in dB/cm/MHz) and thicknesses. 
According to the above equation, we thus assume, as a first approximation, that the attenuation 
increases linearly with the frequency: a(f) = a f. 

The measured backscatter coefficient reported in this study was computed as 

BSC meas (k) = BSC ref (k)^^. (3) 

Pref(k) 

In Eq. (3), the mean backscattered power spectrum P meas was obtained by averaging the power 
spectra of 400 backscattered echoes from blood. The mean power spectrum P re f was obtained from 
a reference sample of non-aggregated RBCs at a low hematocrit of 6% (i.e. Rayleigh scatterers). 9 
In this case, 400 echoes were also averaged. The backscatter coefficient of this reference sample 
BSC re f was estimated using the "Rayleigh estimation" approach used by Yu and Cloutier 1 , which 
theoretical value is given by Eq. (13) in Ref. 8 (three dimensional Perkus Yevick packing factor 
for cylinders). This reference sample was used to compensate the backscattered power spectrum 
Pmeas for the electromechanical system response, and the depth-dependent diffraction and focusing 
effects caused by the US beam. 
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The packing factor W, aggregate diameter D and total attenuation along the propagation path 
a were determined by matching the measured BSC meas given by Eq. © with the theoretical 
BSCtheor given by Eq. ([I]). For this purpose, we searched values of W, D and «o minimizing 
the cost function F(W, D, a )=\ \BSC meas — BSC t heor\ | 2 - hi all studied cases, the cost function 
seemed to have a unique global minimum, as was observed by plotting the cost function surface 
F(W, D) with varying values of olq. An example is given in Fig. 1. 

3. Methods 

3.1 Experimental set up 

US measurements were performed in a Couette flow system to produce a linear blood velocity 
gradient at a given shear rate (see figure 1 in Ref. 11). The system consisted of a rotating inner 
cylinder with a diameter of 160 mm surrounded by a fixed concentric cylinder of diameter 164 
mm. A 60 mL blood sample was sheared in the 2 mm annular space between the two coaxial 
cylinders. The US scanner (Vevo 770, Visualsonics, Toronto, Canada) equipped with the RMV 
710 probe was used in M-mode for this study. The single-element focused circular transducer had 
a center frequency of 25 MHz, a diameter of 7.1 mm and a focal depth of 15 mm. We operated 
at a sampling frequency of 250 MHz with 8 bits resolution (Gagescope, model 8500CS, Montreal, 
Canada). The probe was mounted in the side wall of the fixed outer cylinder and was positioned 
to have its focal zone at the center of both cylinders. To ensure ultrasonic coupling, the hole 
within the outer stationary cylinder (containing the probe) was filled with a liquid agar gel based 
mixture. When solidified, this gel was cut to match the curvature of the cylinder to avoid any 
flow disturbance. The gel was a mixture of distilled water, 3% (w/w) agar powder (A9799, Sigma 
Chemical, Saint-Louis, MO), 8% (w/w) glycerol and a specific concentration of 50 /im cellulose 
scattering particles (S5504 Sigmacell, Sigma Chemical, Saint-Louis, MO) that determined the 
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attenuation coefficient. Four experiments were performed with four mixtures having Sigmacell 
(SC) concentrations varying from 0% to 0.75% (w/w). The 0% concentration constituted the non- 
attenuating gel and the three other mixtures mimicked skin attenuations. 

3.2 Attenuation measurements 

The attenuation coefficients of the reference (0%) and three skin-mimicking phantoms a p were 
determined by using a standard substitution method. Two transducers with center frequencies of 25 
MHz (Vevo 770, Visualsonics, Toronto, Canada) and 20 MHz (V317-SM Panametrics, Waltham, 
MA) were aligned coaxially facing each other for transmission measurements. Transmitted signals 
were recorded both with and without the agar gel sample in the acoustic path. The attenuation 
coefficient was then estimated using a log spectral difference technique. 10 For a given concentra- 
tion of SC, measurements were obtained from two different sample thicknesses, and for each, four 
regions were scanned for averaging purpose. Values obtained were 0.007 ± 0.002, 0.101 ± 0.028, 
0.208 ± 0.029 and 0.317 ± 0.039 dB/cm/MHz for SC concentrations of 0, 0.25, 0.50 and 0.75%, 
respectively. The thickness of the skin-mimicking phantoms e p being fixed to 1 cm, their attenua- 
tion coefficients were thus in the same range as the human dermis (0.21 dB/MHz at 14 - 50 MHz 
considering a 1 mm dermis thickness 12 ). 

3.3 Blood preparation and measurement protocol 

Fresh porcine whole blood was obtained from a local slaughter house, centrifuged and the 

plasma and buffy coat were removed. Two blood samples were then prepared: (i) a H6 reference 

sample, which was a 6% hematocrit non-aggregating RBCs resuspended in saline solution; and (ii) 

a 40% hematocrit T40 test sample, which consisted of RBCs resuspended in plasma to promote 

aggregation. The H6 sample was sheared at 50 s -1 and coupled with the 0% SC concentration agar 

gel. Echoes were selected with a rectangular window of length 0.8 mm at four depths every 0.2 mm 
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(i.e. with 75% overlap between windows). For each depth, the power spectra of the backscattered 
echoes were averaged over 400 acquisitions to provide P re f. Then, the H6 sample was removed 
and the T40 blood was introduced in the Couette device. In the first 30 s, a shear rate of 500 s _1 
was applied to disrupt RBC aggregates. The shear rate was then reduced to residual values of 2, 
5, 10, 20 and 30 s" 1 for 90 s. After that, for each shear rate, acquisitions of 400 rf lines were 
performed for 80 s. Echoes were windowed as for the H6 sample at the same depths and their 
power spectra were averaged to obtain P m eas- This protocol was repeated four times with the four 
agar-based phantoms. 

3.4 Reference measurements with the 0% SC concentration phantom 

The experiment with the 0% SC phantom was realized in order to have reference results on 
packing factors W re f and aggregate diameters D re f obtained from the classical SFSE. 1 These 
parameters were assumed to be true values of packing factors and aggregate diameters for all shear 
rates, and will be compared in the next section with packing factors and diameters estimated by 
the SFSAE and by the SFSE when skin-mimicking phantoms are used. 

It is important to emphasize the fact that the H6 reference sample was also measured with 
the 0% SC phantom. The phantom attenuation, although small with no SC, therefore affected 
equivalently both spectra P meas and P re f in Eq. (3). The resulting measured backscatter coefficient 
BSC re f was thus not biased by attenuation. The terminology "no attenuation" was used for this 
experiment in the following. 

4. Results and discussion 

Figure 2a reports results on W re f and D re f for the SFSE in the case of no attenuation. Typical 
results of the SFSAE minimization procedure for the different agar phantoms at a shear rate of 5 
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s _1 are given in Fig. 2b. All results on W, D and a from the SFSAE are summarized in Fig. 
3 for all residual shear rates. In this figure, the relative errors for each parameter correspond to: 

(W — W re f) /W re f, (D — D re f)/D re f and (a — a re f)/a re f, with a re f measured in transmissions. 
More specifically, a re f corresponds to ^ct^ = ( a p e p + a Uood^biood), where a p e p is the skin- 

i 

mimicking phantom attenuation estimated in transmission, and auood^Uood is the blood attenuation 
taken equal to 0.022 dB/MHz 1 for all shear rates. To underline the necessity to take into account the 
attenuation, parameters W nocomp and D nocomp were evaluated with the SFSE without attenuation- 
compensation when skin-mimicking phantoms were used. Because of the frequency-dependent 
distortion produced by the attenuating medium, large relative errors can be seen in Fig. 4a for both 
parameters. However, by compensating the backscatter coefficients in the SFSE with the value 
measured in transmission (section 3.2), relative errors in Fig. 4b are largely reduced to values 
comparable to those estimated with the SFSAE (see Fig. 3b). 

The SFSAE (Fig. 3) gave quantitatively satisfactory estimates of W, D and a with relative 
errors below 22%, for shear rates between 5 and 20 s" 1 . The SFSE with attenuation-compensation 
(Fig. 4b) gave estimates of W comp and D comp with relative errors below 12% for shear rates between 
2 and 10 s" 1 , and below 28% for the shear rate of 20 s _1 . However, for the SFSAE, the average 
estimates for the shear rate of 2 s _1 were less accurate (relative errors below 57% for W and below 
30% for ao). The estimation of D was satisfactory at that shear rate (relative errors below 14%). 
The worse results of W, D and a were obtained at 30 s -1 for the highest attenuation. 

The apparent limit of applicability of the SFSAE method for shear rates of 2 and 30 s^ 1 may 
be explained by considering the following. At 2 s _1 , for the frequency bandwidth considered (9 
- 30 MHz), the SFSE and consequently the SFSAE seem to reach their limit of applicability for 
large aggregate sizes (typically D re f = 17.5 in Fig. 2a, i.e. kR = 4.8). This limit is illustrated 
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by the bad fit of the SFSE model in Fig. 2a at 2 s -1 . The bad estimations of the SFSAE at 30 
s _1 are explained by the fact that the aggregate diameters were estimated to zero and attenuations 
were overestimated. At this high shear rate, RBC aggregation is partially inhibited and the signal- 
to-noise ratio (SNR) of our measurements was reduced -4 dB between 20 and 30 s _1 for all 
phantoms). The accuracy of the estimates was thus degraded with increasing attenuations, as can 
be seen from the large relative errors at the highest attenuation with the SFSAE but also with the 
SFSE with attenuation-compensation (W comp and D comp ). 

To conclude, the SFSAE performed well for kR < 2.7 (i.e. D = 10 at 5 s" 1 ) and under the 
condition that the SNR is sufficiently good (SNR > 28 dB corresponding to the SNR at 30 s -1 for 
the 0.25% SC). Although the SFSAE gave less accurate estimates for 2 and 30 s -1 , the estimated 
parameter values presented in Fig. 3a show that the SFSAE gave qualitatively satisfactory esti- 
mates for the three SC skin-mimicking phantoms at all shear rates, since the estimates of W and 
D versus shear rates had the same behaviors as W re f and D re f. 

5. Conclusions 

The performance of the new SFSAE was assessed with experimental measurements on blood in 
a Couette flow device. The accuracy of the estimates obtained with the SFSAE was not as satisfac- 
tory as those obtained with the SFSE with attenuation-compensation (i.e when a priori are known 
about the attenuation). Nevertheless, the SFSAE has the major advantage to be easily applicable in 
vivo because of the simultaneous estimation of the blood structural properties and total attenuation 
(contrary to the SFSE attenuation-compensation method, needing the attenuation and thickness of 
the tissue intervening layers to be known). This work thus confirms the in vivo applicability of 
RBC aggregate size and structure estimations. Complementary studies are nevertheless required 

to determine the validity domain of the SFSAE according to kR and attenuation. 
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Figure captions 

Figure 1 . (Color online) (a) Typical aspect of the logarithm of the cost function F(W, D, a ) 
for a fixed value of a®. The logarithm is shown here in order to enhance the visual contrast. 
This cost function has one minimum denoted (W*,D*) that depends on a . (b) Typical 
aspect of the function log (F(W*, D*, a )) for varying values of a (W* and D* being cal- 
culated for each a ). This cost function has a single minimum. 

Figure 2. (Color online) (a) Backscatter coefficients for blood sheared at different resid- 
ual shear rates and measured with the 0% SC concentration phantom (no attenuation), and 
corresponding fitting with the classical SFSE with no compensation for attenuation, (b) 
Backscatter coefficients for blood sheared at 5 s" 1 and measured with each of the four phan- 
toms. The corresponding fitted models are the SFSE for the 0% SC phantom, and the SFSAE 
for the three other skin-mimicking phantoms (0.25, 0.5 and 0.75% SC). 

Figure 3. (Color online) (a) Values of W, D and a (in dB/MHz) for different residual shear 
rates estimated by the classical SFSE for the 0% SC concentration and by the SFSAE for the 
three skin-mimicking phantoms, (b) Corresponding relative errors. 

Figure 4. (Color online) Relative errors of the packing factor and aggregate diameter for the 
three skin-mimicking phantoms obtained with the SFSE (a) with no compensation for atten- 
uation (W nocomp and D nocomp ), and (b) with attenuation-compensation using the attenuation 
values estimated in transmission (W cornp and D cornp ). Parameters W nocomp and W cornp and 
similarly D nocomp and D comp are compared with W re f and D re f, respectively. 
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